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Recent studies have shown the possibility of water transport across carbon nanotubes, even in the case of 
nanotubes with small diameter (0.822 nm). In this case, water shows subcontinuum transport following an 
ordered ID structure stabilized by hydrogen bonds. In this work, we report MD simulations describing the 
effect of a perpendicular electric field in this single-file water transport in carbon nanotubes. We show that 
water permeation is substantially reduced for field intensities of 2-3 V/nm and it is no longer possible under 
perpendicular fields of 4 V/nm. 

1 Introduction 

Despite its strong hydrophobic character, single walled carbon nanotubes (SWCNTs) are spon- 
taneously filled by water [1, .2|. Transport of water in these systems has novel and promising 
properties which have been extensively studied by Molecular Dynamics (MD) simulations (see 
for example [3] for a review) . Quite interestingly, the water transport trough a carbon nanotube 
depends strongly on its radius and MD simulations predict a transition between continuum to 
subcontinuum transport as the radius decreases [1], as recently confirmed experimentally [H] . 
In the case of the nanotubes with small radius, the nanotubes are filled by a one-dimensional 
ordered chain of water molecules which maintain 2 hydrogen bonds per molecule inside the CNT 
[6]. This effective ID water system has extremely interesting properties and it has been studied 
experimentally [5], by MD simulations [THU EHH] and analytical treatments based on the ID 
Ising model [9]. Interestingly, this system has been analyzed as a model for water transport in 
biological pores fTOj. 

A recent MD simulation study has shown the strong response of this single-file water system 
inside CNTs to electric fields parallel to the axis of the nanotube [8], due to the peculiarities of 
the ID hydrogen bonding network. Motivated by this previous study, in this work we consider a 
related but different question, namely the effect of applying an electric field perpendicular to the 
single file water transport across carbon nanotubes. In this case, a competition arises between 
the hydrogen bonding network (which tends to maintain the alignment of the chain of water 
molecules in the nanotube axis direction) and the perpendicular electric field (which tends to 
align the water dipole in the perpendicular direction). A simple numerical calculation proposed 
in [11] suggests that this effect will appear for field intensities between 2-3 V/nm. Noting that a 
single water- water hydrogen bond has a free energy of about bksT and the dipole of the water 
molecule in liquid phase is about 2.5 D we can say that H-bonding will compete with external 
fields with intensities up to ~ 2.4 V/nm. Albeit high, these field intensities can be studied 
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in simulations and experiments. In this work, we report a preliminary set of MD simulations in 
order to study the effect on water transport of this competition between hydrogen bonding and 
electric field. 



2 Methods 

We have performed a total of 33 molecular dynamics (MD) simulations of water transport across 
carbon nanotubes in both equilibrium and nonequilibrium conditions in presence and in absence 
of an electric field perpendicular to the axis of the nanotubes. 

All simulations described here were performed using the 2.7 version of the NAMD program |12j . 
The preparation of the system and analysis of the results were made using the Visual Molecular 
Dynamics (VMD) software |13j . The simulated system is shown in Figure [TJ It consisted of a 




Figure 1. Representative snapshot of MD simulations. Atoms from water molecules are shown with Van der Waals radius 
and carbon atoms from the four SWCNT are shown with lines. Note that water molecules fill the tubes as single file with 
a strong orientation. 

periodic simulation box of dimensions 23.0 x 19.9 x 30.4 containing 248 water molecules and a 
porous membrane made of four carbon nanotubes with their axis aligned in the z direction. Each 
carbon nanotubes has 144 carbon atoms and has a radius of 0.411 nm and length 1.34 nm. Water 
molecules were modelled using the TIP3P model as implemented in the CHARMM force field and 
carbon atoms from nanotubes were modelled as Lennard-Jones spheres with parameters e = 0.07 
kcal/mol and a = 3.9848A . Carbon atoms were maintained fixed during all the simulation as in 
previous works[5l[6]. Lennard-Jones interactions were computed using a smooth (10-12 A) cutoff, 
as it is customary done in NAMD2 simulations. The electrostatic interactions were calculated 
using the particle-mesh Ewald (PME) method. We employ a multiple time step. The equations 
of motion are solved with a 1 fs time step, nonbonded interactions are updated each 2 time steps 
and electrostatic interactions are updated each 4 time steps. The temperature was maintained 
at 300 K in all simulations employing a Langevin thermostat with a 5 ps~^ dumping constant. 
The Langevin thermostat is applied to the water oxygen atoms only. 

Some of our simulations were performed in presence of an external pressure difference. This 
external pressure was applied in NAMD2 as a Tel force (as implemented in |14j) which acts on 
oxygen atoms in a slab of 5.4 A thickness (all oxygen atoms with coordinates z > 12.5 A or 
z < —12.5 A). This methodology is also described in detail in Ref |15| . The applied pressure AP 
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is computed from the force / by AP = nf /A where n is the number of water molecules in the 
slab and A its transversal area. 

In most of our simulation runs, we applied an uniform electric field Ey, perpendicular to the 
axis of the nanotubes, a feature which is implemented in NAMD2. The field is acting in all the 
simulation box. 

In all simulations, we compute the number of permeation events, defined as the number of 
water molecules observed to completely cross the nanotube (entering from one end and leaving 
the tube from the other end). Obviously this quantity depends on the simulation time, so it was 
normalized to the total production time (in ns). The length of the production runs was decided 
by monitoring this quantity, i.e. simulations were performed until this quantity reached an stable 
average value. Depending on the particular simulation runs, simulation times are between 5 and 
40 ns. 



3 Results and Discussion 



First, we conducted equilibrium simulations (no applied external fields or pressure gradients). 
In this case, the average flux of water is negligible, as expected in equilibrium. In equilibrium 
we obtain an average rate of permeation events of water molecules of 9 molecules/ns in any of 
the directions. 
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Figure 2. Permeation events of water molecules per ns observed in MD simulations as a function of an applied electrical field 
perpendicular to the axis of the nanotubes. Open triangles correspond to equilibrium simulations and permeation events 
are averaged over the +2 and —2 direction. Filled triangles correspond to simulations with a water flow induced under an 
external pressure of AP = 3714 atm. In this case, permeation events occur only in the direction of decreasing pressure. 



Now, we consider the effect of applying an external electric field Ey perpendicular to the axis 
of the nanotubes. The results are shown in Figure [2] Fields up to 1 V/nm have no significant 
impact on water transport. However, for larger fields, the number of permeation events decreases 
substantially as the field strength increases. For fields about 2 V/nm, the number of permeation 
events has been reduced to less than its equilibrium value (see Fig. [2]). We also observe that 
chains inside nanotubes start to break and we even observe crossing of isolated (not hydrogen 
bonded) water molecules across the nanotube, as shown in the snapshot in Figure [3j Larger 
fields reduce the permeation events even more and it can be say that water permeation is not 
possible for fields larger than 4-5 V/nm. 

We have also performed simulations under a large external pressure AP = 3714 atm (see Figj2] 
). In absence of external electrical field, we observe about 133 permeation events per ns. In all 
cases permeation occurs in the direction of decreasing pressure, as should be expected in this 
case. Under the application of external, perpendicular electric field we obtain results strongly 
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Figure 3. Snapshot of MD simulations under an electric field of 2 V/nm in the y direction. Atoms from water molecules 
are shown with Van der Waals radius and carbon atoms from the four SWCNT are not shown for simplicity. Note that one 
of the tubes is empty and in another nanotube there is only a single, isolated water molecule inside one of the nanotubes 
oriented in the field direction. In another of the tubes the hydrogen-bonded water chain is still intact and water molecules 
fill the tube. 



resembling to those obtained in the absence of a pressure gradient. Again, the effect of fields 
smaller than 1 V/nm is negligible, and substantial reduction of the flow of water requires fields 
of the order of 2-3 V/nm. It can be said that the presence of the external pressure has no effect 
on the competition between hydrogen bonding and the external electric field. 



4 Conclusions 



In this work, we have performed MD simulations of single-file transport of water in carbon 
nanotubes of small radius. Simple arguments suggest that the application of an electric field 
perpendicular to the axis of the carbon nanotube will disrupt the hydrogen bond structure of 
this ID water system and affect the transport properties of water inside the nanotube. Our 
simulations show that water transport is not affected under perpendicular fields up to 1 V/nm. 
Larger fields substantially reduce the flow of water, and in fact transport is not possible for 
perpendicular fields above 4 V/nm. Our work shows the possibility of using perpendicular elec- 
tric fields as nanofluidic "switches" to regulate water transport in carbon nanotubes. It is also 
interesting to note that a previous work [8j has shown the advantages of pumping water across 
nanotubes using parallel electric fields (between 0.01 V/nm - 1 V/nm) instead of using a pres- 
sure gradient. It could be interesting to analyze also this parallel case with higher fields (2-4 
V/nm) which could also somewhat influence the hydrogen bonding network. Our own ongoing 
work indicates that under parallel fields ~ 4 V/nm there are significant structural changes in the 
single file of water molecules. Overall, all these results show the interesting possibilities arising 
from using electric fields to control water flow in carbon nanotubes. 



June 23, 2011 



0:16 Molecular Simulation nanotubes'23May 



Molecular Simulation 5 



Acknowledgments 

This work is supported by the Spanish Government (grants FIS2009-13370-C02-02 
and CONSOLIDER-NANOSELECT-CSD2007-00041) and GeneraUtat de Catalunya (grant 
2009SGR164). 



References 

[1] K. Wu, B. Zhou, P. Xiu, W. Qi, R. Wan and H. Fang, Kinetics of water filling the hydrophobic channels of narrow 
carbon nanotubes studied by molecular dynamics simulations, J. Chem. Phys. 133 (2010), art. num. 204702. 

[2] H. Kumar, B. Mukherjee, S. Lin, C. Dasgupta, A. K. Sood and P.K. Maiti, Thermodynamics of water entry in 
hydrophobic channels of carbon nanotubes, J. Chem. Phys. 134 (2011), art. num. 124105. 

[3] A. Alexiadis and S. Kassinos, Molecular Simulation of Water in Carbon Nanotubes, Chem. Rev. 108 (2008), pp. 
50145034 

[4] J. A. Thomas and A. J. H. McGaughey , Water Flow in Carbon Nanotubes: Transition to Subcontinuum Transport, 

Phys. Rev. Lett. 102 (2009), art. num. 184502 
[5] X. Qin, Q. Yuan, Y. Zhao, S. Xie and Z. Liu, Measurement of the Rate of Water Translocation through Carbon 

Nanotubes, 11 (2011) 2173-2177 
[6] G. Hummer, J. C. Rasaiah and J. P. Noworyta, Water conduction through the hydrophobic channel of a carbon 

nanotube. Nature 414 (2001), pp. 188-190. 
[7] G. Zuo, R. Shen, S. Ma and W. Guo, Transport properties of single-file water molecules inside a carbon Nanotube 

biomimicking water channel, ACS Nano 4 (1010), pp 205-210. 
[8] J. Su and H. Guo, Control of Unidirectional Transport of Single-File Water Molecules through Carbon Nanotubes in 

an Electric Field, ACSnano 5 (2011) 351-359. 
[9] J. Kofinger and C. Dellago, Single-file water as one- dimensional Ising model. New J. Phys. 12 (2010) art. num. 093044. 
[10] F. Zhu and K. Schulten Water and Proton Conduction through Carbon Nanotubes as Models for Biological Channels, 

Biophys. J. 85 (2003) pp 236244. 
[11] J. Faraudo The missing link between the Hydration Force and interfacial water: Evidence from computer simulations. 

Current Opinion Colloid and Interf. Sci. doi:10.1016/j.cocis.2011.04.009 (accepted for publication) 
[12] J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, and K. 

Schulten, Scalable molecular dynamics with NAMD, J. Comp. Chem. 26 (2005), 1781-1802. 
[13] W. Humphrey, A. Dalke, and D. J. Schulten, VMD - Visual Molecular Dynamics, Molec. Graphics 14 (1996), 33-38. 
[14] J. Cohen, F. Zhu and E. Tajkhorshid, Simulation of Water Permeation Through Nanotubes, NAMD tutorial available 

at http:/ /www. ks.uiuc.edu/Training/Tutorials/ 
[15] F. Zhu, E. Tajkhorshid and K. Schulten, Pressure- Induced Water Transport in Membrane Channels Studied by Molec- 
ular Dynamics, Biophys. J. 83 (2002) 154 160. 



